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Abstract 

Background: Upon co-stimulation with CD3/CD28 antibodies, activated CD4 + T cells were found to lose their 
susceptibility to HIV-1 infection, exhibiting an induced resistant phenotype. This rather unexpected phenomenon 
has been repeatedly confirmed but the underlying cell and molecular mechanisms are still unknown. 

Methods: We first replicated the reported system using the specified Dynal beads with PHA/IL-2-stimulated and 
un-stimulated cells as controls. Genome-wide expression and analysis were then performed by using Agilent whole 
genome microarrays and established bioinformatics tools. 

Results: We showed that following CD3/CD28 co-stimulation, a homogeneous population emerged with uniform 
expression of activation markers CD25 and CD69 as well as a memory marker CD45RO at high levels. These cells 
differentially expressed 7,824 genes when compared with the controls on microarrays. Series-Cluster analysis 
identified 6 distinct expression profiles containing 1,345 genes as the representative signatures in the permissive 
and resistant cells. Of them, 245 (101 potentially permissive and 144 potentially resistant) were significant in gene 
ontology categories related to immune response, cell adhesion and metabolism. Co-expression networks analysis 
identified 137 "key regulatory" genes (84 potentially permissive and 53 potentially resistant), holding hub positions 
in the gene interactions. By mapping these genes on KEGG pathways, the predominance of actin cytoskeleton 
functions, proteasomes, and cell cycle arrest in induced resistance emerged. We also revealed an entire set of 
previously unreported novel genes for further mining and functional validation. 

Conclusions: This initial microarray study will stimulate renewed interest in exploring this system and open new 
avenues for research into HIV-1 susceptibility and its reversal in target cells, serving as a foundation for the 
development of novel therapeutic and clinical treatments. 
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Background 

A comprehensive picture of the host factors putatively 
supporting HIV-1 (human immunodeficiency virus, type 
1) replication in cells has emerged from recent siRNA 
studies [1-4] and meta-analysis [5]. Genome-wide land- 
scapes of host genes and proteins involved in HIV-1 infec- 
tion and disease progression have also been established in 
gene array studies [6-9] and novel proteomic approaches 
[10-12]. In parallel, genome-wide association studies 
(GWAS) have revealed a set of inheritable genetic varia- 
tions in large populations related to susceptibility to HIV- 
1 infection [13-15]. 

Against these exciting developments, there is still no 
global view of the host cellular factors that render the 
target cells resistance to infection, in spite of a few well- 
studied restriction factors [16-19]. A recent genome- 
wide screening for novel restriction factors [20] further 
highlighted this awareness and our interest in establishing 
a more holistic pictures of the host determinants working 
against HIV-1 susceptibility. 

A crucial starting point for studies of this type is the 
consideration of sample sources, cell types, and experi- 
mental settings. In this regard, CD4 + T cells are the first 
choice; they are the major cell type amongst all the sus- 
ceptible targets and reservoirs of HIV-1 infection [21]. 
Apart from their intrinsic susceptibility, CD4 + T cells' ac- 
tivation in vivo (during the natural courses of HIV-1 infec- 
tion) and in vitro (typically with PHA/IL-2 stimulation) is 
generally recognized as an absolute prerequisite for the 
virus to replicate productively [22]. However, rather unex- 
pectedly, Levine et al. (1996) found that activation by co- 
stimulation with CD28 led to a complete loss of suscepti- 
bility to HIV-1 infection in these cells. This phenomenon 
was subsequently confirmed by several independent 
groups [23-26]. In spite of the fact that several further 
studies attempted to explore the potential of using thus- 
stimulated cells for the treatment of SIV/HIV infection in 
monkey models and clinical settings, progress in elucidat- 
ing the underlying cellular and molecular mechanisms 
seems to have halted since 2002 for unknown reasons. 

In the present study, in order to investigate the pos- 
sible mechanisms of the observed reversal of HIV-1 
susceptibility in these activated CD4 + T cells, we first 
replicated the reported experimental settings [27] and 
then performed genome-wide expression analysis using 
Agilent micr oar rays. 

Methods 

Isolation and stimulation of CD4 + T cells 

Buffy coats were obtained from healthy donors (Kunming 
Blood Station) and peripheral blood mononuclear cells 
(PBMCs) were isolated by Ficool-Hypaque (TBD Sciences) 
gradient centrifugation. Resting CD4 + T cells were then 
purified by magnetic negative selection (CD4 + T cell 



Isolation Kit II, Miltenyi Biotec) and only the aliquots of 
cells with purity > 95% as determined by flow cytometry 
were used for further analysis. This study was reviewed 
and approved by the internal review board of the Kunming 
Institute of Zoology, Chinese Academy of Sciences (ap- 
proval ID: RTYX20090910-1, approval date: 2009-09-10). 
All donors provided written informed consent for partici- 
pation in this study. 

Cells were stimulated as previously reported [27,28]. 
Briefly, freshly isolated CD4 + T cells were resuspended 
in RPMI 1640 medium (Gibco) supplemented with 10% 
heat inactivated fetal bovine serum (Gibco) and 20 mM 
HEPES (Amresco) and seeded in 6-well plates either at 
an initial density of 2 x 10 6 cells/ well with human recom- 
binant IL-2 (100 U/ml) and PHA (5 ug/ml) or 0.5 x 10 6 
cells/well with polystyrene beads coated anti-CD3/CD28 
antibodies (Dynal beads CD3/CD28 T Cell Expander, 
Dynal) at a bead to cell ratio of 3:1. Half media were 
changed every 2 days in CD3/CD28 costimualted CD4 + T 
cells and every 3 days in PHA/IL-2 stimulated cells. Cells 
were then cultured at 37°C in a humidified incubator with 
5% C0 2 for 6 days. Cells and the derived RNA samples 
and data sets were labeled "P", "R", and "B", respectively, 
according to the PHA/IL-2 stimulated, un-stimulated rest- 
ing and beads-stimulated settings. 

Flow cytometry 

Aliquots of cells to be analyzed were washed with FACS 
buffer (PBS supplemented with 1% bovine serum albumin 
(BSA)) and stained for 40 minutes in the dark at 4°C 
with anti-CD45RO-FITC, anti-CD25-PE, anti-CXCR4, 
fluorescent dye 5-(and -6)-carboxyfluorescein diacetate 
succinimidyl ester (CFSE) (Sigma-Aldrich), anti-Ki67 
(Abeam), anti-CCR5 (Biolegend) and anti-CD69-PE-Cy- 
5 (BD Biosciences) and appropriate isotype controls. 
After washing with FACS (Fluorescence Activated Cell 
Sorter) buffer, cells were fixed with 4% paraformalde- 
hyde and analyzed on a FACS Calibur using Cell Quest 
and Flowjo 7.6.1. 

Gene expression profiling with Agilent microarrays 

Total RNA was extracted from stimulated or un-stimulated 
CD4 + T cells by TRIzol (Invitrogen) followed by a purifica- 
tion using RNeasy columns (Qiagen) according to the man- 
ufacturers protocols. The amount and quality of RNA 
preparations were evaluated on an Agilent 2100 Bioanalyzer 
with RNA6000 Nano Reagents and Supplies (Agilent). 

Quality-checked RNAs were then transcribed with the 
First-Strand cDNA Synthesis Kit (Agilent) and their 
expression data obtained using Agilent 4 x 44 K Human 
Whole-Genome 60-mer oligonucleotide microarrays 
according to the protocols by the manufacturer. The ori- 
ginal microarray data from this study are available at the 
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NCBI GEO database (http://www.ncbi.nlm.nih.gov/geo/) 
under the accession number GSE34252. 

Analysis of microarray data 

The normalized ratio of the gene expression signals was 
log 2 transformed and hierarchical clustering was per- 
formed with average linkage. The clustered heatmap 
was visualized using Treeview. The RVM (Random vari- 
ance model) f-test was applied to filter differentially 
expressed genes for the different situations. After the 
significance analysis and FDR (false discovery rate) 
analysis, we selected the differentially expressed genes 
according to the p-value and FDR threshold set at p < 
0.01 and FDR <0.01 [29-31]. And the fold changes of 
any two groups are more than 2. 

Series-Cluster analysis was performed to identify the 
global trends and model profiles of expression according 
to signal density under the "P", "R" and "B" conditions 
and in the P-R-B sequence. Fishers exact test and the 
multiple comparison test [32,33] were applied to identify 
the model profiles with probability significantly higher 
than expected as random. 

Gene ontology (GO) analysis [34] was performed to fa- 
cilitate elucidating the biological implications of unique 
genes in the significant or representative profiles. GO 
analysis was used to find the main function of the genes 
having the same expression trend according to the Gene 
Ontology. Fisher s exact test and ^ test were applied to 
identify the significant GO categories and FDR was used 
to correct the ^-values. 

Gene co-expression networks analysis [35] were per- 
formed to track the interactions among the differentially 
expressed genes, according to their normalized signal in- 
tensity in 6 representative profiles. Pearson correlation 
was applied to each pair of genes and the significantly 
correlated pairs were used to construct the network [36]. 
To locate the core regulatory genes in the networks, k- 
core scoring was introduced to simplify graph topology 
analysis [37,38]. A k-core of a given gene indicates its 
hub or nodal status with connection to "k" other genes 
in a network [37,38]. Accordingly, the genes with largest 
k-core scores and highest degrees of connection were 
identified as "key regulatory genes" in a network [39] 
and those unique to each network were selected as 
"marker genes" for the cellular status of intrinsic suscep- 
tibility or induced resistance to HIV-1 infection. 

Results 

Cellular characterization of stimulated and un-stimulated 
CD4+ T cells 

Following the previously published details of the initial 
experiment (Levine et al. 1996) and other independent 
studies [23-26], we re-established de novo and verified 
the experimental systems (Figure 1). After stimulation 



for 6 days with the CD3/CD28 antibodies coated beads 
(labeled "B" for beads), purified CD4 + T cells became en- 
larged and highly proliferative, forming a large amount of 
big cell colonies (Figure 1A). When compared with the 
un-stimulated (labeled "R" for resting) or PHA/IL-2-stim- 
ulated (labeled "P" for PHA) cells at the functional level, 
the induced resistance to HIV-1 infection in these co- 
stimulated cells [23-28] was also readily replicated in our 
earlier experiments (described as a brief report in a home 
journal [40]). We monitored the cell viability and prolifer- 
ation first by Trypan blue exclusion test which showed 
more than 95% viable cells for all the experiments (data 
not shown) and then by Ki67 and CFSE staining 
(Additional file 1: Figure SI and Additional file 2: Figure 
S2). We then proceeded to analyze the surface pheno types 
of the generated cells and found that a homogeneous 
population of cells emerged following the co-stimulation 
(Figure IB). In contrast to the controls, these cells were 
uniformly and simultaneously stained positive for both the 
general and early activation markers (CD25 and CD69, re- 
spectively) at high levels (Figure IB). Another important 
marker for activated and memory T cells, CD45RO, was 
also uniformly expressed following co-stimulation in these 
"B" cells but in a much lower percentage than in "R" or 
"P" cells (Figure IB). 

Global expression profiles and signatures 

Using Agilent 4 x 44 K Human Whole-Genome Micro- 
arrays, we found that the overall gene expression pat- 
terns were clearly and sharply different among the "R", 
"P" and "B" cells. As shown in Figure 2 (left) with their 
clustered heatmaps and hierarchical patterns, in total 
7,824 genes represented 8,128 transcripts (out of the 
41,000 test probes on the chips) were differentially 
expressed according to the RVM (Random variance 
model) algorithm (p-value < 0.05, FDR < 0.05) (Additional 
file 3: Table SI). 

In subsequent Series -Cluster analysis, we identified 16 
possible profiles (Figure 2 middle) (Additional file 4: 
Table S2), which represent the overall expression pat- 
terns. Of these, 10 (#1, 5, 6, 7, 8, 9, 10, 11, 12 and 15, 
Figure 2 middle) showed the sharp difference between "B" 
and "P" cells. More specifically, genes down-regulated in 
"P" cells are shown in the "Pdown" (#1) profile and those 
up-regulated in these cells were grouped in the "P up " (#9) 
profile (Figure 2 right). Similarly, genes up- or down- 
regulated only in "B" are represented in "B up " (#8) and 
"Bdown" (#15) profiles. We merged three profiles (#5, 6 
and 7) that contain the genes down-regulated in "P" but 
up-regulated in "B" into a new "PdownB up " profile, in order 
to show their special importance in the reduced resistance 
as they contain the genes which were down regulated in 
the susceptible "P" cells and at the same up-regulated in 
the resistant "B" cells. In the same way, we merged profiles 
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Figure 1 Morphology and surface phenotypes of stimulated and un-stimulated CD4+ T cells. A: CD4 + T cells stimulated with PHA/IL-2 
("P") or CD3/CD28 coated beads ("B") were observed on day 0 (un-stimulated ("R")), day 3 and day 6. B: Flow cytometric analysis of cell surface 
phenotypes. "R", "P" and "B" cells were stained with antibodies to CD25 (FITC), CD69 (PE-Cy5) and CD45RO (FITC). Numbers indicate the 
percentage of each subset. A majority population of cells emerged in "B" cells which expressed CD25, CD69 and CD45RO. We did the 
experiments in 19 biological replicates. 



#10, 11 and 12 into a new "P up Bdown" profile to show the 
genes that were up-regulated in "P" cells but at the same 
time down-regulated in the "B" cells (Figure 2 right). 

These 6 representative profiles contain 1,345 differentially 
expressed genes in "B" and "P" cells (Figure 3 left) and give 
the apparent expression signatures either unique to "P" and 
"B" cells or polarized expressed in "P" and "B" cells at two 
opposite extremities. The remaining 6,479 genes differen- 
tially expressed in the other 6 profiles (#2, 3, 4, 13, 14 and 
16, Figure 2 middle) were unique to the "R" cells, up- or 
down-regulated when compared with "P" and "B" cells as 
the baseline. Due to the nature of the "R" cells that are 
seemingly susceptive to the infection but do not support 
active viral replication, the biological importance of these 



genes in HIV-1 infection are more complicated to dissect 
and interpret. We kept aside these 6,479 "R" cell-unique 
genes for future studies. 

Functional categories and significant ontologies of the 
differentially expressed genes in "P" and "B" cells 

It is immediately noticeable (Figure 3 left) that most of the 
1,345 differentially expressed genes in the 6 representative 
profiles falls in the "B" cell-related expression, with 471 
genes in the "B down " profile and 465 in "B up ", respectively, 
representing the down-regulated and up-regulated un- 
permissive genes in the resistant "B" cells (Figure 3 left, 
Additional file 5: Table S3). In the "B down " profile, 222 genes 
with known functions gave a rather complicated ontology 
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Figure 2 Overall patterns of 7,824 differentially expressed genes and the 6 representative profiles in the permissive "P" and resistant 
"B" cells. Patters were plotted on the heatmap using Treeview. Red represents up-regulated genes while green represents down-regulated 
genes. Hierarchical clustering is shown on the left. All 16 expression profiles identified by Series-Cluster analysis are shown in the middle and 6 
summarized representative profiles are shown on the right. 



category but, after applying Fisher s exact test, j 2 test, and 
FDR analysis, a clearer picture emerged. Statistically, 43 
genes were significantly involved in 4 major functional 
categories (Figure 3 middle), showing: (a) response to 
stimulus (immune response, cell surface receptor linked 
signal transduction and positive regulation of natural killer 
cell mediated cytotoxicity); (b) metabolic process; (c) 
multicellular organismal development; and (d) cellular 



process are predominantly influenced by the co- 
stimulation. 

Other part of exactly the same 4 major categories was 
found significantly influenced in the "P up " and "P up Bdown" 
profiles, with part of one more category (localization) being 
also up-regulated in the permissive "P" cells or down- 
regulated in the un-permissive "B" cells (Figure 3 middle). 
Combined GO analysis of the genes in the "P up ", "P up Bdown" 
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Figure 3 Gene ontology (GO) analysis and significant functional genes. Left: Numbers of genes in the 6 representative profiles. Middle: 
significant functional categories and the number of genes in them. Each color represents one category and the size of each sector in a pie 
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and B down " profiles showed one more category (cell prolif- 
eration) was also significantly influenced (Figure 3 right). 
Altogether, these 6 significantly affected functional categor- 
ies contain 101 genes (Figure 3 right) representing the 
whole set of the "permissive genes" that are associated posi- 
tively with the highly susceptible status or negatively with 
the induced resistance of the cells. As shown in Figure 4A, 
these permissive genes are clearly hierarchical and highly 
enriched in the categories related to defense, immune/in- 
flammatory response and signal transduction under the 
general heading of response to stimuli [41]. 

More genes (352) with known functions were identified 
in the "B up " profile with 117 being significantly involved 
in 7 major GO categories, showing that (a) cellular process 
(cell adhesion, cell differentiation and regulation of tran- 
scription), (b) multicellular organismal development, (c) 
metabolic process, (d) response to stimulus, (e) localization, 
(f ) circulatory system process, and (g) others are predomin- 
antly up-regulated following the co-stimulations. Exactly 
the same 7 categories were also up-regulated in the 
"P down B up " profile, although the numbers of the genes in 
each category varied (Figure 3 middle right). These 7 cat- 
egories contain all the 4 significant ones in the "P down " 



profile and one more category (cell fate determination) was 
also identified. Altogether, the 8 categories contain 144 
genes representing the whole set of the "resistant genes" 
that are associated positively with the induced resistance or 
negatively with the highly susceptible status of the cells 
(Figure 3 right). However, they are dispersed in more 
diverse categories (Figure 4B) than the permissive ones 
(Figure 4A), indicating that the resistant genes have wider 
influence in cellular functions and may hold the crucial 
check points for HIV-1 infection in "multicellular organis- 
mal development", "cellular process" and "metabolic 
process". One of the striking examples is that 30 genes 
involved in the process of cell adhesion were found up- 
regulated in "B" cells (profiles "B up " and "PdownB up ") 
(Figure 4B) (Additional file 6: Table S4). 

Direct interactions of the genes with known functions 
with HIV-1 proteins 

When subjected to further screening, 22 of the 144 po- 
tential "resistant genes" and 23 of the 101 potential 
"permissive genes" were identified for their direct inter- 
actions with HIV-1 proteins (Figure 5 A & Additional file 7: 
Table S5), according to the data in the "HIV-1 Human 




Figure 4 Hierarchical GO categories of genes involved in induced permissive "P" and resistant "B" cells. GO categories trees were 
hierarchically built using the Gene Ontology Enrichment Analysis Software Toolkit (GOEAST): http://omicslab.genetics.ac.cn/GOEAST/). The category tree 



of response to stimulus related genes in "P up ", "B down " and "P up B down " cells is shown in A. Category trees of multicellular organismal development, cellular 

B up " and "P down B up " cells are shown in B. 



process and metabolic process related genes in "P d( 
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Figure 5 (See legend on next page.) 
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(See figure on previous page.) 

Figure 5 Interactions of significant functional genes with HIV-1 proteins. Filled circles represent genes involved in significant functional 
categories. Red represents genes in profiles "P d0 w n ", " B u P " and "P d0 wnB U p", while green represents genes in profiles "P up ", "B d0 w n " and "P up B d0 wn". 
Unfilled squares represent HIV-1 proteins. A: Overall interactions of significant functional genes with HIV-1 proteins; B: Functional categories of 
known gene interaction with HIV-1. 



Protein Interaction Network" (http://www.ncbi.nlm.nih. 
gov/RefSeq/HIVInteractions/). Their significant functional 
categories are shown in Figure 5B. 

Specifically, 3 (Env, Nef, and Tat) of HIV-1 proteins 
attracted most of the host (no matter permissive or resist- 
ant) factors in this set of data, suggesting that host factors 
may have the capacity to counterbalance with the permis- 
sive genes by interacting with the same early genes in HIV- 
1 infection. This is coincided with the initial report that the 
novel anti-HIV effect induced by co-stimulation happens 
around viral entry and before integration [27]. However, in 
terms of host gene functions, no overlap was identified be- 
tween the resistant and permissive genes as the former 
were mapped on the categories of cell adhesion, multicellu- 
lar organismal development and cell differentiation, while 
the latter on immune response, cell surface receptor signal 
transduction, cellular defense response and chemotaxis. 
Fewer but prominent genes, e.g., CD59, IL-2, IL-15, INDO 
and FN1 were found to interact with Gag and Vpr, al- 
though we did not find any genes in all our profiles capable 
of interacting with HIV-1 Vif or Rev. 

Co-expressed genes and their networks in the permissive 
"P" and un-permissive "B" cells 

Given that HIV-1 tends to interact with "key" host pro- 
teins, such as bottlenecks and hubs in gene interaction 
networks [42,43], we used the genes in the representa- 
tive profiles and constructed 6 gene co-expression net- 
works (Figure 6 left). We then applied the "k-core" 
scores (see Methods) to identify those with highest net- 
working degrees as the "key regulatory" genes that may 
play pivotal roles in gene interactions and regulations. 
Thus, 53 genes in profiles "P up ", "B down " and "P U pB d own"> 
which tentatively contain the "permissive genes" were 
chosen as "key regulatory" genes and, among them, 29 
genes with known functions (underlined in Figure 6 
middle) are involved in immune response, signal trans- 
duction and so on (Figure 6 right). Similarly, 84 genes in 
profiles "Pdown"> "B up " and "PdownB up "> which tentatively 
contain the "resistant genes", stood out with the highest 
k-core scores and, out of them, 36 genes (underlined in 
Figure 6 middle) are involved in cell adhesion, multicel- 
lular organismal development and so on (Figure 6 right). 

KEGG pathways and schematic overview of HIV-host gene 
interactions in the resistant "B" cells 

We then mapped all the genes with known functions in 
the significant categories (Figure 3) and all the "key 



regulatory" genes (Figure 6) on KEGG pathways from the 
Kyoto Encyclopedia of Genes and Genomes database. A 
clear picture emerged in the resistant "B" cells showing 
(Figure 7) that (a) the up-regulated genes in "B up " and 
"PdownB U p" profiles were enriched in the cellular processes 
(Figure 8) including: filamentous actin (F-actin), tight 
junction, actomyosin assembly contraction, proteasomes, 
proteolysis, lysosomes, degradation, and Na + Ca 2+ ex- 
change, and (b) the down-regulated genes in "B down " and 
"PupBdown" profiles were enriched in the process of actin 
polymeration, apoptosis, cell cycle checkpoint, ER to Golgi 
transport and terminally misfolded. 

Verification of the differential gene expression by FACS 

We then performed FACS analysis to examine the expres- 
sion levels of two HIV-1 co-receptors, CCR5 and CXCR4. 
As shown in (Figure 9), CCR5 was up-regulated in the per- 
missive "P" cells but down-regulated in the resistant "B" 
cells, while CXCR4 remained down-regulated in both "P" 
and "B" cells. These data not only validated our microarray 
analysis — albeit only a small part of the whole data set — but 
also were consistent with the initial findings that the in- 
duced resistance was only effective to CCR5-dependent 
strains of the HIV-1 virus [24,44,45]. This also confirmed 
our previous report [40] showing that, following co- 
stimulation, the levels of CCR5 expression varied to some 
degrees in the 12 individuals tested, but the general trend 
was clearly down-regulated. Conceivably, the down- 
regulation of the critical co-receptor CCR5 accounts, at least 
in part, for the induced resistance in the co-stimulated cells. 

Clues for novels genes accounting for HIV-1 susceptibility 
and resistance 

It should be noted that the aforementioned co-expression 
network analysis was performed using all the genes with 
and without known functions. Strikingly, among the total 
137 "key regulatory" genes, only 9 genes have been previ- 
ously reported to interact with HIV-1 (Table 1) while the 
remaining majority (128 genes) as of yet have no confirm- 
ation of their potentials to interact with HIV-1 genes or pro- 
teins. When going back to the 1,345 differentially expressed 
genes in the 6 representative profiles (Figure 2 left), around 
a third of them (513 genes) do not yet have any assigned 
functions. In view of their prominence in the representative 
profiles and co-expression networks, these unknown genes 
constitute a new source for future studies on their roles in 
biological processes and potential involvement in deter- 
mining HIV-1 susceptibility and resistance. 
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Figure 6 Co-expressed genes and their networks. Left: Number of genes in 6 representative profiles (Figure 3 middle left). Red nodes 
represent "key regulatory" genes while blue nodes represent other regulated genes. Node size represents the power of the interrelation among 
the nodes, and edges between two nodes represent interactions between genes (i.e. the more edges of a gene, the more genes connecting to 
it, the more central role it has within the network). Middle right: Name list of "key regulatory" genes with highest k-core. Green represents 
unknown genes; known functional genes were underlined. Right: significant function of these genes. 
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Discussion 

In the present study, we re-established an experimental 
system for cellular and molecular studies on the induced 
resistance to HIV-1 infection in CD4 + T cells as first 
reported by Levine et al [27]. This cellular system includes 
two methods most commonly used to activate and expand 
CD4 + T cells in vitro. Incubation with PHA/IL-2 induces 



proliferating CD4 + T cells which are highly susceptible to 
HIV-1 infection and highly permissive for the subsequent 
viral replication, whereas co- stimulation with CD3/CD28 
reverses the intrinsic susceptibility in these cells and ren- 
ders them with resistance to HIV-1 infection and un- 
permissiveness for viral replication. Thus, these two 
extremely polarized statuses of activated CD4 + T cells, 
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Figure 9 Verification of CCR5 and CXCR4 expression by flow cytometry. Cells of "R", "P" and "B" were stained with antibodies to CCR5 (PE) 
and CXCR4 (PE) as described in Methods. Numbers indicate the percentage of gated subset that expressed CCR5 or CXCR4. These results have 
been replicated and confirmed in 10 individuals. 
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Table 1 Interactions between HIV-1 proteins and 9 "key 
regulatory" genes 

Profiles Genes HIV-1 Mode of interactions 



proteins 



B UP FN1 


Env 


binds 




Gag-Pol 


cleaves 




Tat 


competes with, modulated by, 
up-regulates 




Nef 


up-regulates 


LAMB3 


Tat 


up-regulates 


MYH10 


Gag-Pol 


cleaves 


TGM5 


Env 


interacts with, modified by 


VEGFA 


Env 


modulates 




Tat 


cooperates with, induces release of 


Bdown PRKCZ 


Env 


interacts with, inhibits, regulated by, 
up-regulates 




Gag-Pol 


inhibited by, phosphorylated by 




Tat 


Activates, regulated by, 
phosphorylated by 


PdownBup FLJ22814 


Nef 


up-regulates 


PupBdown MGAT5B 


Env 


processed by 


SMPD3 


Env 


activates 



with the un-stimulated resting cells (susceptible to the in- 
fection but not permissive for the rapid viral replication) 
as control in the middle of the "spectrum of permissive- 
ness", constitute a highly insightful system in the search 
for host factors responsible for HIV-1 susceptibility and 
its reversal in CD4 + T cells. 

We were attracted more by the biological significance 
and potentials of the induced resistance in the "B" cells. 
We first replicated this in vitro system and verified its use- 
fulness (Figure 1 and [40]) at the cellular level. Considering 
the controversy around the source cells, i.e., which subsets 
(naive or memory CD4 + T cells) could be induced to gen- 
erate the resistance [24,26,46], we adopted a strategy using 
total CD4 + T cell population to start our study without 

Table 2 Genes involved in cell cycle arrest 



Profiles Genes of cell cycle arrest GenBank accession 

B up APBB2 NM_1 73075 

B up CDKN1A NM_078467 

B up & PdownBup DST NM_01 5548 

B up GAS2L3 NM_1 74942 

B up NAPSA NM_004851 

B up NOL3 NM_003946 

B up SESN2 NM_031459 

B up TCF4 NM_003199 



prior cell sorting before co-stimulations. This gave us the 
chance to address the issue of which subset of CD4 + T 
cells was generated to be responsible for the induced resist- 
ance and led us to the first set of our findings, which shows 
a homogenous cell population emerged from the co- 
stimulation. These cells expressed three activation markers 
(CD25, CD69 and CD45RO) at high levels at the same 
time. Although the expression of these classical activation 
markers may not ultimately account for the induced resist- 
ance, they served well as indicators for the highly activated 
population and, more relevant to the following analysis, 
the homogeneity of the induced cells gave us the confi- 
dence for the subsequent whole-genome wide search for 
the "key regulatory genes" and the "core marker genes" in 
the status switch to resistance. 

Apart from the above novel findings in defining the cell 
population emerged during the reversal of susceptibility, 
the present study is the first whole-genome-wide analysis 
for the genes that account for the induced resistance. Once 
again, the system allowed us to pin down the potential per- 
missive and resistant genes since susceptible "P" and the re- 
sistant "B" cells are polarized at the two extremities of the 
"permissive spectrum". As a result, the 1,345 genes in the 6 
specific profiles for the "B" and "P" cells could be truly the 
representatives accounting for the intrinsic susceptibility to 
HIV-1 infection and its induced reversal in the target cells. 
This simple dichotomy classification further allowed us to 
analyze the differentially expressed genes for their gene 
ontology and co-expression with ease and confidence. 

We have not included the popular pathway and tran- 
script factor analysis in this initial genome-wide search 
but, by focusing on the gene categories and co-expression, 
managed to identify the major cellular processes involved 
the reversed susceptibility in the major HIV-1 target cells. 
Prominently, these include actin cytoskeleton system, pro- 
tein degradation and cell cycle arrest. First, several lines of 
evidence highlighted that actin cytoskeleton were regarded 
as a barrier and the hijack of actin cytoskeleton facilitates 
entry of HIV into its target cells [47,48]. In our study, 



Description 

Homo sapiens amyloid beta (A4) precursor protein-binding, family B, 
member 2 (Fe65-like) (APBB2), mRNA [NM_1 73075] 

Homo sapiens cyclin-dependent kinase inhibitor 1A (p21, Cipl) (CDKN1A), 
transcript variant 2, mRNA [NM_078467] 

Homo sapiens dystonin (DST), transcript variant 1eA, mRNA [NM_0 15548] 

Homo sapiens growth arrest-specific 2 like 3 (GAS2L3), mRNA [NM_1 74942] 

Homo sapiens napsin A aspartic peptidase (NAPSA), mRNA [NM_004851] 

Homo sapiens nucleolar protein 3 (apoptosis repressor with CARD domain) 
(NOL3), mRNA [NM_003946] 

Homo sapiens sestrin 2 (SESN2), mRNA [NM_031459] 

Homo sapiens transcription factor 4 (TCF4), mRNA [NM_003199] 
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8 genes (ACTN2, PTPRF, SORBS1, PAR3, PVRL2, MYL9, 
MYH10 and PAK2) were found up-regulated in the resist- 
ant "B" cells ("B up " profile in Figure 8A). These genes me- 
diate assembly and contraction of actin and actomyosin 
and formation of tight junctions. Whether they formed a 
barrier against HIV-1 infection warrant further functional 
verifications. To the contrary, genes associated with actin 
polymerization (FN1 and ITGAX, Figure 8A), which have 
been shown to promote viral binding and entry [48] and is 
also necessary in chemotaxis and cytokinesis [49] was 
sharply down regulated in the "B" cells. Given that chemo- 
taxis may actually serve to fuel the infection response by 
recruiting susceptible, activated CD4 + T cells to the virus, 
ultimately aiding viral dissemination [49], we therefore hy- 
pothesized that co-stimulation could also block the tracks 
for HIV-1 release and dissemination by reducing actin 
polymerization and chemotactic response. 

Second, genes involved in ubiquitin-dependent degrad- 
ation by the proteasome (WNT3, FZD1, FZD2, FZD3, 
FBXW11, RAD23B, DBAHV2, HLPA4L, IRAK3, IL1RAP, 
PROC and PROS1), which is the third step of endoplasmic 
reticulum associated protein degradation, were also pre- 
dominantly in "B" cells (Figure 8B). This is consistent with 
previously report that the ubiquitin-proteasome was nega- 
tively associated with HIV-1 replication by acting to des- 
troy incoming viral complexes at the early steps [50,51]. 
We also found that ATPVOA4 and IGF2R, which are 
important for lysosome-mediated degradation of endo- 
cytosed proteins [52,53], were up-regulated in "B" cells 
(Figure 8B). Therefore, co- stimulation could enhance the 
cellular defense to HIV-1 infection by up-regulating these 
lysosome- and proteasome-related genes [54]. 

Third, it is also noticeable that 8 unique genes (Table 2) 
involved in cell cycle arrest were up-regulated in "B" 
cells. Apparently, these genes may have helped for the 
co-stimulated cells to reach the highly proliferative sta- 
tus. However, their role in HIV-1 /AIDS infection and 
disease progression need further detailed studies since 
modifications to the cell cycles was observed in both 
long-term non-progressors (LTNPs) and progressive pa- 
tients [55,56]. Anyhow, this line of evidence reflects that 
host target cells may have intrinsic mechanisms, which, 
when activated, would halt viral infection at the specific 
point of the cell cycle. 

The other two major cellular programs which were 
profoundly influenced by the co-stimulation were me- 
tabolism and apoptosis (Figure 4B). This seems incon- 
sistent with previous studies showing that (a) numerous 
metabolism-associated genes were down-regulated in 
LTNPs [55] and (b) apoptosis is widely accepted as a 
mechanism for T cell depletion in vitro and in vivo 
[3,57,58]. Further 5 genes (CA2, MED28, PLOD2, SLC2A3 
and TMTC1), which were up-regulated as probably 
resistant genes in our "B" cells, were proposed as "host 



dependent factors (HDF)" in recent small-interfering 
RNA-knockdown screens [1-4,20]. More surprisingly per- 
haps, in our set of potentially resistant genes, we did not 
find any overlap with those reported to be associated with 
LTNPs [55], viral latency [59,60], or even the "resistance 
genes" identified by virus-host interaction network ana- 
lysis using public data sets [61] (data not shown). How- 
ever, these apparent discrepancies do not necessarily 
discredit our studies but, to the contrary, may favorably 
underscore the uniqueness and importance of the highly 
controlled experimental system we adopted and call for 
further studies down the line. 

Finally, we chose to examine the expression of CCR5 
and CXCR4, the 2 major HIV-1 co-receptors, on the cell 
surface by FACS analysis as the first part of our ongoing 
efforts to verify our array data and analysis. To our satis- 
faction, their expression in this set of experiments and 
our previous screening [40] consistently corroborated 
our bioinformatics analysis and helped in explaining the 
fact that the induced resistance by the co-stimulation is 
a CCR5-dependent phenomenon [24,44,45]. Obviously, 
large-scale validation and functional studies at both 
mRNA and protein levels are needed to corroborate the 
role of the key regulatory genes (e.g., ITGAX, VEGFA, 
FN1, CCND1, CA12, LAMB3, MYL9, etc.) and their path- 
ways as the major determinants for the induced resistance. 
We hope the present study serves as a fresh call for a 
renewed interest in the experimental system as a useful 
model for the search of novel host resistant factors. 

Conclusions 

In summary, we replicated a simple yet powerful cellular 
system with CD3/CD28 co-stimulation and confirmed 
its usefulness in studying the induced resistance to HIV- 
1 in CD4 + T cells. This initial microarray study, al- 
though still descriptive and correlative in nature, allowed 
us a chance to glean valuable new insights into this 
phenomenon. Based on the overall expression patterns 
and apparent signatures of the differentially expressed 
genes, we managed to pin down 245 (101 potentially 
permissive and 144 potentially resistant) significant func- 
tional genes and 137 (84 potentially permissive and 53 
potentially resistant) "key regulatory" genes involved in 
the reversal of target cell susceptibility to HIV-1 infec- 
tion. Moreover, we showed that this system has the 
potential as a rich source for the search of novel genes 
accounting for the intrinsic susceptibility to HIV-1 and 
its reversal in the major target cells. We hope our find- 
ings will stimulate renewed interest in investigating the 
mechanisms underlying this phenomenon. We surmise 
that future studies of this kind will help in furthering 
our understanding of HIV-1 infection, which may ultim- 
ately then lead to the development of novel biomarkers 
and therapeutics. 
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Additional files 



Additional file 1: Figure SI. Comparative CFSE staining among "R", "P" 
and "B" cells. Cells of "R" were stained with CFSE and cells of "P" and "B" 
were stained with CFSE on day 3 and day 6. Significant T cell 
proliferation was observed during the culture period. 

Additional file 2: Figure S2. Comparative Ki67 staining among "R", "P" 
and "B" cells. Cells of "R" were stained with anti-Ki67 and cells of "P" and 
"B" were stained with anti-Ki67 on day 3 and day 6 and it also showed 
significant! cell proliferation during the culture period. Red line: isotype 
control; green line: Ki67. 

Additional file 3: Table SI. Differentially expressed genes. 

Additional file 4: Table S2. Profiles of differentially expressed genes. 

Additional file 5: Table S3. Gene ontologies of the differentially 
expressed genes in "P" and "B" cells. 

Additional file 6: Table S4. Hierarchical GO categories of genes 
involved in induced permissive "P" and resistant "B" cells. 

Additional file 7: Table S5. Interactions between HIV-1 proteins and 
significant functional genes. 
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